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STATISTICAL  METHODS  \?PR  OPR  I ATE  PCR  EVALUATION  OP  FUZE 
EXPLOSIVE -TRAIN  SAPErj.  aRD  RELIABILITY 

Prepared  by: 

H.  P.  Culling 


ABSTRACT:  The  statistical  problems  of  evaluating  fuze  explo- 
sive trains  for  safety  and  reliability  are  presented.  The 
Bruceton  and  Probit  procedures  of  testing  are  illustrated  by 
examples  of  actual  tests  and  the  interpretation  of  the  results 
Is  discussed. 

Appendix  A presents  a step-by-step  description  of  the 
calculations  necessary  to  determine  the  "50$  reliable  distance", 
the  standard  deviation,  errors  of  estimate,  and  "degree  of 
reliability"  by  the  Bruceton  Method  of  Analysis.  Tables  and 
graphs  necessary  for  the  calculations  are  appended. 

Appendix  B presents  the  Probit  Method  calculations  neces- 
sary to  determine  the  Probit  Line  Equation  from  which  estimates 
of  the  "50$  unsafe  distance",  standard  deviation,  errors  of 
estimate,  and  "degree  of  safety"  may  be  obtained.  A test  for 
"goodness  of  fit"  of  the  Probit  Line  Is  presented.  Tables 
necessary  for  the  analysis  are  appended. 

Appendices  A and  B thus  provide  a working  handbook  for 
the  two  methods  of  analysis,  which  are  bot*  appropriate  for 
either  safety  or  reliability  tests. 

The  methods  presented  and  tables  appended  have  been 
extracted  from  references  (a)  and  (b)  and  adapted  to  the 
problems  of  estimating  fuze  safety  and  reliability.  Physical 
considerations  governing  the  applicability  of  the  test  methods 
described  to  various  configurations  of  explosive  trains,  and 
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the  use  of  the  results  In  the  evaluation  of  ordnance  Items 
as  a whole,  arc-  the  subject  of  a separate  report  and  are 
not  considered  herein. 

Many  similar  problems  in  ordnance  evaluation  and  testing 
in  which  the  data  must  be  of  a "go-no  go"  nature  may  pro- 
fitably be  investigated  by  the  methods  presented  in  this  report. 
For  example,  the  minimum  and  average  firing  energy  charac- 
teristics of  electric  primers  are  often  so  determined. 

Since  the  test  points  are  distributed  over  the  center  of 
the  distribution  rather  than  being  concentrated  at  the  nd 
points  (where  all  or  none  of  a small  sample  may  be  expected 
to  function  and  which  usually  correspond  to  nominal  part- 
dimensions,  firing  voltages  or  whatever),  these  methods  are 
relatively  efficient  in  affording  a good  basis  for  estimating 
the  effects  of  later  design  or  manufacturing  changes  or 
proposals.  The  value  of  the  data  from  tests  conducted  in 
accordance  with  these  plans  is  therefore  less  likely  to  be 
lost  in  case  of  such  changes. 

The  applicability  of  these  methods  is  dependent  upon  an 
understanding,  based  upon  experience  with  large  numbers  of 
previous  tests  or  demonstration  of  the  validity  on  a theoret- 
ical basis,  of  the  nature  of  the  physical  processes  involved. 

Xn  particular,  since  these  methods  are  aimed  at  giving  an 
overall  estimate  of  the  population  from  a very  limited  number 
of  samples,  the  analysis  cannot,  in  itself,  check  the  degree 
to  which  the  assumptions  of  normal  distribution  and  linearity 
(or  other  assumed  relationship)  between  the  two  variables  are 
justified.  The  greatest  value  of  the  test  methods  is  in  - 
quickly  confirming  the  suitability  of  a design  or  pointing 
out  a marginal  condition  requiring  correction  or  further 
investigation  by  the  use  of  more  samples. 
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13  October  1953 


This  report  has  been  prepared  under  Task  N0L-B2b-4l-l , for 
the  following  purposes: 

1.  To  illustrate  the  statistical  philosophy  of  evaluating 
fuze-explosive  trains  for  safety  and  reliability, 

2.  To  present  the  application  of  two  methods  of  analysis 
which  deal  with  "go-no  go"  results  as  a function  of  a con- 
trolled variable,  and 

3»  To  provide  a wo.  ing  handbook  for  using  these  methods  in 
the  design,  analysis  and  presentation  of  explosive-train 
safety  and  reliability  evaluations. 

These  procedures  have  been  developed  by  the  Trains  and  Mecha- 
nisms Branch,  Mechanical  Evaluation  Division,  in  the  course 
of  the  evaluation  of  ordna.nce  items  under  various  Bureau  of 
Ordnance  tasks.  Tables  appearing  in  this  report  have  been 
taken  or  adapted  froi..  the  references  cited. 

Tile  Naval  Ordnance  Laboratory  is  indebted  to  Professor  Sir 
Ronald  A.  Fisher,  Cambridge,  to  Dr.  Frank  Yates,  Rothamsted, 
and  to  Messrs.  Oliver  and  Boyd  Ltd.,  Edinburgh,  for  per- 
mission to  reprint  Table  No.  IX  from  their  book  "Statistical 
Tables  for  Biological,  Agricultural,  and  Medical  Research". 


JOHN  T.  HAYWARD 
Captain,  USN 
Commander 
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STATISTICAL  METHODS  APPROPRIATE  FOR  EVALUATION  OF  FUZE 
EXPLOSIVE -TRAIN  SAFETY  AND  RELIABILITY 


INTRODUCTION 

1.  Present-day  emphasis  on  highly  reliable  and  safe  weapons 
calls  for  Increasingly  accurate  methods  of  estimating  the  pro- 
babilities of  reliable  and  safe  functioning  of  components. 

The  usefulness  of  statistical  procedures  in  determining 
reliability  is  discussed  in  references  (f)  and  (g) . The 
general  procedure,  however,  is  based  on  a large  sample  of 
quantitative  data  which  is  obtained  for  a given  test  con- 
dition. This  report  illustrates  the  statistical  methods 
which  are  applicable  to  "go-no  go"  data  (such  as  the  firing 
or  non-firing  of  an  explosive  element)  obtained  at  various 
levels  of  a test  condition. 

2.  Organizations  engaged  in  the  evaluation  of  fuzes  to 
determine  their  suitability  for  release  to  production  have 
the  problem  of  determining  the  reliability  and  safety  of  fuze 
explosive  trains.  This  report,  which  is  limited  to  the 
statistical  aspects  of  the  problem,  has  the  following 
objectives: 

a.  To  present  the  statistical  philosophy  of  evaluating 
fuze  explosive  trains  for  safety  and  reliability, 

b.  To  illustrate  the  application  of  two  methods  of 
analysis  (the  Bruceton  and  Probit)  which  deal  with  "go-no  go" 
results  as  a function  of  a controlled  variable,  and 

c.  To  provide  a working  handbook  for  using  these  methods 
in  the  design,  analysis  and  presentation  of  explosive-train 
safety  and  reliability  evaluations. 


Inadequacy  of  Confidence  Limit  Estimates 

3.  In  order  to  meet  the  explosive-safety  requirement,  an 
explosive  train  usually  contains  an  out-of-line  safety  device 
which  rotates  or  slides  from  the  unarmed  to  the  armed  position. 
The  safety  or  reliability  of  the  train  may  be  estimated  roughly 
by  testing  In  the  unarmed  and  armed  positions,  respectively. 
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But  If  the  sample  size  Is  limited  to,  say,  20  for  each  test, 
the  statistical  interpretation  of  the  results  cannot  establish 
the  high  assurance  of  safety  or  reliability  necessary  for 
release.  For  example,  given  a sample  of  20,  all  of  which 
fire  reliably  in  the  armed  position,  one  may  expect  from 
future  samples  a reliability  as  low  as  83$  (or  17$  failures). 

If  the  sample  is  increased  to  100,  one  may  still  expect  a3 
many  as  4$  failures.  These  estimates  of  percentage  failures 
are  based  on  the  lower  95$  confidence  limit  (reference  (a)). 
Since  an  explosive  train  is  probably  more  reliable  than 
Indicated  by  such  a test,  a more  accurate  estimate  is  necessary. 


Method  of  0 ,ing  Between  Arr  * and  Unarmed  Positions 

4.  If  the  safety  device  were  placed  half  way  between  the 
armed  and  unarmed  positions  and  the  explosive  train  still 
functioned  reliably,  one  would  expect  the  reliability  of  the 
train  in  armed  position  to  be  considerably  better  than  shown 
in  the  previous  example.  Then,  if  more  samples  were  tested 
at  various  distances  of  partial  arming  and  in  such  a manner 
that  an  array  of  high  to  low  percentage-functioning  distances 
was  obtained,  a graph  could  be  plotted  of  percentage  reliable 
versus  distance  from  armed  position.  By  methods  of  statis- 
tical extrapolation,  a percent-reliable  curve  could  be  calcu- 
lated from  which  an  estimate  of  the  99$-rellable  distance 
could  be  read.  If  the  99$~reliable  distance  occurs  when  the 
safety  device  is  still  considerably  out  of  line,  it  Is 
probable  that  the  explosive  train  will  function  consistently 
In  the  armed  position. 

5.  The  safety  test  could  be  run  In  a similar  fashion,  where 
an  array  of  high  to  low  percent-unsafe  conditions  versus 
distance  was  obtained.  The  calculated  percent-unsafe  curve 
would  give  the  estimated  l$-unsaf*e  distance. 


Methods  of  Analysis  and  Calculated  Parameters 

6.  Since  the  estimates  of  the  99$-reliable  distance  and 
l$-unsafe  distance  are  the  bases  of  conclusions  concerning 
reliability  and  safety  of  the  train,  the  method  of  extra- 
polation for  obtaining  these  distances  should  be  as  accurate 
as  possible.  The  method  should  also  be  applicable  for  small 
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samples.  These  two  restrictions  prohibit  the  use  of  several 
established  methods  which  employ  relatively  simple  test  pro- 
cedures. Two  established  methods  which  meet  these  require- 
ments are  the  Bruceton  Method  of  Sensitivity  Testing  (reference 
(b))  and  Probit  Analysis  (reference  (c)).  Both  methods  may 
be  employed  to  estimate  the  50^-reliable  distance,  Xg,  or  the 
50^-u.nsafe  distance,  Xy,  and  the  respective  standard  deviations, 
op  and  on.  Each  X minus  its  appropriate  a gives  an  estimate 
of  the  84^  distance.  If  0 is  added  to  X,  an  estimate  of  the 
15$  distance  is  obtained.  By  using  the  appropriate  constants 
obtained  from  a table  of  normal  curve  areas,  one  may  estimate 
any  percentage  distance.  Therefore,  the  estimates  of  XR  and 
5Ty  and  or  and  Cy  are  the  only  measurements  necessary  to 
determine  a measure  of  reliability  or  safety.  In  appendices 
A and  B,  detailed  descriptions  of  the  calculation  of  these 
measurements  for  the  Bruceton  Method  and  the  Probit  Method 
are  presented. 


Specific  Applications  for  Each  Method 

7.  Although  both  of  the  methods  provide  a means  of  estimating 
the  same  measures  (XR,  Xy,  aR,  ay),  there  are  specific  reasons 
for  the  application  of  one  instead  of  the  other,  depending  on 
the  purpose  of  the  test. 

a.  The  Bruceton  Method  of  testing  centers  the  testing  In 
the  region  of  the  50^  distance , while  the  Probit  Method  of 
testing  includes  the  high  and  low  percentage  points.  Thus, 
the  Bruceton  Method  Is  more  accurate  for  estimates  of  XR  and 
Xy,  while  the  Probit  Method  is  more  accurate  for  estimates  of 
oR  and  Oy. 

b^  The  Bruceton  Method  requires  fewer  samples  to  deter- 
mine Xp  or  Xy  and  a rough  estimate  of  oR  and  ay,  so  it  may  be 
used  for  small  samples  (not  less  than  20) . It  should  be  held 
in  mind  that  the  Bruceton  estimates  of  aR  and  ay  tend  to  be 
slightly  less  than  the  Probit  estimates;  and,  in  connection 
with  safety  or  reliability  tests,  such  estimates  are  overly 
optimistic . 

c.  The  Bruceton  Analysis  may  not  be  used  except  when  the 
Bruceton  scheme  has  been  followed,  using  equally  spaced  test 
distances  and  the  up-and-down  procedure  of  testing  (described 
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in  appendix  A) . On  the  other  hand,  the  Probit  Analysis  is 
adaptable  to  data  of  unequal  sample  sizes  at  each  test  distance 
and  the  test  distances  need  not  be  equally  spaced. 

d.  The  Bruceton  Method  of  Analysis  is  based  on  the 
assumption  that  the  percentage  functioning  is  related  to  tv 
test  distances  by  means  of  the  normal  error  function,  while 
the  Probit  method  may  be  analyzed  by  other  appropriate  trans- 
formations, if  necessary. 

Thus,  the  Probit  Method  may  be  considered  more  adaptable  to 
most  test^  data^  except  when  the  sample  sizes  are  too  small, 
or  when  Xr  or  Ay  are  the  most  important  measures.  In  some 
cases  the  Bruceton  Method  of  Analysis  is  preferred  because 
the  computations  Involved  are  relatively  simple,  but  this 
fact  alone  should  not  determine  the  choice  of  method.  See 
reference  (d)  for  a more  detailed  discussion  of  the  two 
methods . 


Comparison  of  Reliability  or  Safety  Probabilities 

8.  The  purpose  of  an  explosive-train  reliability  or  safety 
test  is  to  determine  Just  how  reliable  or  safe  the  train  is 
in  the  armed  or  unarmed  positions,  respectively.  It  Is 
possible  by  means  of  the  proposed  methods  to  determine  the 
probability  of  each;  but,  in  most  cases,  this  probability  is 
in  the  range  of  99$  to  99*999999$.  Obviously,  small  samples 
do  not  warrant  so  many  significant  figures,  yet  the  trains 
represented  by  many  nines  after  the  decimal  point  are  cer- 
tainly more  reliable  than  those  with  only  one  or  two.  Another 
disadvantage  of  the  probability  figure  Is  that  It  doe3  not 
Indicate  the  amount  of  variation  due  to  sample  differences  in 
the  explosive  components  and  mechanical  parts  (or  inconsistency 
in  manufacture  of  the  design).  In  order  to  determine  a 
difference  between  two  high  probabilities,  a measure  which 
gives  the  ratio  of  X to  a should  be  used. 


A Measure  of  the  Degree  of  Reliability  or  Safety 

9.  Two  measures  termed  the  "degree  of  reliability".  Dr,  and 
the  "degree  of  safety",  Dg,  which  are  not  subject  to  tne 
discrepancies  of  the  probability  values,  may  be  calculated. 
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Dr  Is  defined  as  Xr  in  terms  of  or,  or 


where  Xn  Is  measured  from  the  armed  position  for  reliability. 
Ds  is  obtained  in  a similar  manner  by 


where  X™  is  measured  from  unarmed  position.  The  values  for 
X and  0 are  obtained  from  the  Bruceton  or  Probit  estimates  of  . 
the  reliability  curve  and  safety  curve,  respectively.  A com- 
parison of  two  reliabl-lity  functions  which  have  identical  Xr's 
but  different  values  of  cR  may  best  Illustrate  the  usefulness 
of  this  measure . 


The  above  graph  shows  XR  for  I and  II  is  at  2.5/32  an  inch. 
The  reliability  curve  I has  a aR  equal  to  0.5/32  Inches,  and 
curve  II  has  a aR  equal  to  1/32  Inches.  The  reliability  pro- 
bability (P)  in  the  armed  position  for  curve  I equals  99-99997# 
while  curve  II  has  a P = 99.4#.  In  a practical  sense,  these 
two  values  of  P do  not  indicate  any  significant  difference  when 
estimated  from  small  samples.  But  by  means  of  the  degree 
measure,  we  find  that  curve  I has  a Dr  = 5-0  while  curve  II 
gives  a Dr  = 2.5.  By  comparing  the  two  Dr  values,  one  can 
conclude  that  DR  for  the  train  of  curve  I is  twice  the  value 
of  Dr  for  the  train  of  curve  II,  while  the  probability  figures 
might  lead  to  the  conclusion  that  the  two  trains  are  almost 
similar. 
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10.  A more  difficult  comparison  between  probability  values 
occurs  when  two  curves  have  slightly  different  values  of  Xr 
and  or,  but  cne  probability  values  occur  in  the  99  to  100% 
region.  Here  again  the  Dr  of  each  will  indicate  which  is  the 
more  reliable . 


Lower  Limit  for  the  Degree  of  Reliability  or  Safety 

11.  Consideration  must  be  given  to  the  degree  of  reliability 
or  safety  that  is  necessary  to  insure  consistent  operation  of 
the  explosive  train.  A lower  limit  for  DR  and  Ds  has  been 
designated  as  5.0  to  insure  better  than  99.5$  operation.  This 
lower  limit  was  calculated  to  allow  for  the  errors  of  estimates 
from  small  samples  of  not  less  than  20.  If  the  exact  50$  dis- 
tance, Xr,  and  the  exact  standard  deviation,  or,  were  known, 
99*5$  operation  would  be  assured  with  a degree  equal  to  2.6. 

(In  other  words,  99^5$  reliability  is  insured  if  Xp/aE  = 2.6 
where  2.6  Is  the  normal-deviate  constant  corresponding  to 
99-5$-)  But  the  small  sample  estimates  of  X and  o are  subject 
to  sampling  errors,  so  these  must  be  allowed  for  in  addition 

to  the  degree  of  2.6. 

12.  The  error  in  X (sample  estimate)  may  be  calculated  by: 

o^  - a/spn,  where  n Is  the  number  of  results  that 
responded,  but  not  the  total  sample  size,  and  o 
is  the  standard  deviation  estimated  from  the 
sample . 

If  we  consider  samples  of  20  or  more  with  at  least  9 favorable 
results,  then  a*  will  be_less  than  or  equa_l  to  a/ JW  or  o/3. 

The  95$  lower  limit  for  X can  be  found  by  X - 1.86o^,  or 

X - - ^2 0 which  equals  X - 0.6a.  Therefore,  to  insure  that 
3 __ 

the  exact  Xr  has  been  included  in  the  measure  of  D,  an  allowance 
of  0.6a  must  be  subtracted  from  the  sample  X. 

This  95$  lower  limit  for  X may  be  explained  by  the 
following  short  discussion.  Confidence  limits  may 
be  computed  for  a mean  by  the  expressions. 
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* ± tos  or  X ± t 2_  (since  cs  = , 

where  the  value  for  t may  be  found  ir.  table  2, 
appendix  A using  n-  - 1 degrees  of  freedom  and 
that  probability  or  '’risk"  of  excluding  X which 
we  are  willing  to  allow.  Thus,  if  n =*  9 (as  in 
the  example  above)  and_we  wished  to  be  90$  sure 
of  including  the  real  X (i.e.,  we  desired  to 
know  the  90$  confidence  limits),  we  should  enter 
the  t table  with  9 - 1 or  8 degrees  of  freedom 
and  a probability  level  of  100$  - 90$  or  10$. 

This  would  give  a t value  -of  1.86,  and  30  the 
90$  confidence  limits  for  n *-=  9 a re 

X ± 1.86  £Lu  or  X ± 0.6a. 
v9 

The  other "way  of  looking  at  these  90$  confi- 
dence limits  is  that  the  probability  is  10$ 
that  the  real  X has  been  excluded . This  10$ 
probability  is  the  sum  of  the  5$  probability 
that  the  true  X lie3  above  the  greater  90$ 
confidence  limit,  X + 0.6a,  and  the  5$  pro- 
bability that  the  true  X lies  below  the 
smaller  90$  confidence  limit,  X - 0.6a.  In 
other  words,  the  probability  is  95$  that  the 
true  X lies  above  the  smaller  90$  confidence 
limit,  X - 0.6a,  which  is  the  value  given 
above  for  the  95$  lower  limit  for  X. 

13.  The  sampling  error  of  a may  be  determined  by  obtaining 
an  upper  limit  of  a which  will  include  the  true  aE  with  95$ 
confidence.  A relation  between  aE  and  the  sample  a Is  given 
by  the  F-ratio,  defined  as 

F - a|/a2. 

Tables  are  available  in  most  statistical  texts  which  give 
values  of  F for  specific  confidence  levels  and  for  various 
sample  sizes  used  for  estimating  a.  The  95$  upper  limit  of 
a is  calculated  using  the  F-value  given  at  the  5$  level 
under  Ni  - 00  (representing  the  size  of  the  population)  and 
jj2  s 8 (representing  n - 1 which  was  used  to  obtain  a) . This 
value  of  F Is  2.9  which  gives  the  95$  upper  limit  of  a as: 
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- 2.9  or  1.7o  = aE. 

Therefore , the  error  in  a is  allowed  for  by  estimating  oF 
as  I.70. 

14.  was  predion  sly  sta'^d  that  99.5#  reliability  was 
insured  if 


^g/ *“  2.6. 

By  substituting  the  errors  for  the  sample  estimates  in  the 
above  relation,  it  becomes. 


X - 0.6a 
TT7a 


2.6  or  X/a  « 5-0. 


Thus  D,  obtained  from  a sample  of  20  or  more,  must  be  greater 
than  5.0,  if  99.570  operation  is  desired.  (It  should  be  noted 
that  the  samples  tested  must  be  representative  of  the  produc- 
tion lots  for  these  figures  to  be  meaningful.) 


Looking  back  on  this  problem,  it  can  be  stated 
that  if  Dp  (which  is  the  ratio  of  X to  a)  is 
5.0  or  greater,  then  even  if  the  3aaq?le  esti- 
mates of  a and  X were  both  so  overly  optimistic 
that  they  might  be  exceeded  only  5#  of  the  time 
from  a sample  as  small  as  9 the  ratio  Xg/og  (the 
true  X to  the  true  0)  should  be 'equal  to  or 
greater  than  2.6.  It  was  stated  above  that-  a 
value  of  2.6  would  insure  at  least  99*5# 
operation  (or  safety).  The  requirement  implied 
here  is  that  the  ratio  of  the  pessimistic  95# 
limits  for  0 and  X be  greater  than  2.6.  The 
value  of  5.0  was  merely  found  to  be  the  necessary 
lower  limit  to  the  ratio  of  the  estimates  of  X 
to  0 from  a sample  as  small  as  9 to  be  certain 
that  the  2.6  requirement  for  the  ratio  of  the 
pessimistic  limits  for  X and  0 was  met.  It  can 
be  seen  that  while  a D.,  of  5.0  always  Implies 
99.5#  operation,  the  failure  of  a set  of  data 
to  yield'  a Dg  of  5.0  does  not  mean  that  the  2.6 
value  for  the  ratio  of  the  pessimistic  95# 
limits  for  X and  0 (and  hence  99*5#  operation) 
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has  not  been  met,  since  the  sample  size  may 
have  been  larger  than  9 and  this  would  result 
in  the  pessimistic  95#  limits  for  X and  o 
being  closer  to  the  sample  estimates  of  3T 
and  a.  The  ratio  of  these  pessimistic  95# 
limits  for  X and  a should  be  computed  when- 
ever Dr  or  Dg  is  found  to  be  less  than  5»0, 
since  if  this  ratio  iB  equal  to  or  greater 
than  2.6,  the  statement  that  the  item  is  99*5# 
safe  or  operable  may  3till  be  made. 


Degree  of  Reliability  Versus  Degree  of  Safety 


15.  Since  an  explosive  train  must  be  both  reliable  and  safe, 
these  two  opposing  requirements  must  be  balanced.  In  other 
words,  Dr  must  not  exceed  Dg  to  such  an  extent  that  the  safety 
is  endangered,  and  conversely.  In  a theoretically  ideal  train, 
the  difference  between  Dr  and  Dg  would  be  minimized . In 
practice,  the  difference  may  be  large  so  long  as  both  values 
exceed  the  lower  limit  of  5.0.  Tables  1 through  7 present 
actual  test  results  determined  by  the  Bruceton  and  Probit 
Methods.  Graphs  1 through  7 are  the  respective  graphs  of 
percent-reliable  and  percent-unsafe  conditions  estimated  from 
the  data  given  in  the  tables.  The  calculated  DR  and  Dg  for 
each  of  the  samples  are  as  follows:  1 0 


% 

Degree  "of 
Reliability 


DS 

Degree" of 
Safety 


Table 

1, 

Graph 

1 

ft 

2, 

IT 

2 

ft 

3, 

If 

3 

n 

4, 

ff 

4 

it 

5, 

ft 

5 

it 

6, 

ft 

6 

it 

7, 

ft 

7 

10.4 

6.8 

3.4* 

12.2 

8.5 

6.6 

3.6* 


14.7 

2.0* 


12.3 

12.0 

23.1 

31.8 


The  tabulation  indicates  that  the  device  of  Graph  2 was 
unsafe  and  those  of  Graphs  3 and  7 were  unreliable;  The 
other  devices  may  be  compared  by  their  values  of  Dr  and  Ds 
which  show  that  the  devices  for  which  Graphs  1 and  4 were 
plotted  contain  the  best  combination  of  safety  and  reliability. 


3 


* 
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CONCLUSIONS 

16  Fuze  explosive-train  reliability  or  safety  estimates 
obtained  by  means  of  confidence  limits  cannot  establish  the 
high  assurance  necessary  for  release.  Therefore,  tests  con- 
ducted with  the  explosive-safety  device  in  the  armed  or 
unarmed  position  are,  for  the  most  part,  inadequate  for  deter- 
mining statistical  estimates  of  reliability  or  safety, 
respectively. 

17.  The  Bruceton  or  Probit  methods  of  testing  at  partial- 
arming  positions  provide  results  from  which  a statistical 
estimate  of  a high  assurance  of  reliability  or  safety  can  be 
established.  If  the  degree  of  reliability  (the  ratio  of  the 
50^-reliable  distance  to  the  standard  deviation)  is  greater 
than  5.0,  the  explosive-train  reliability  will  be  at  least 
99.5^.  Similarly,  if  the  degree  of  safety  (the  ratio  of  the 
50^-unsafe  distance  to  the  standard  deviation)  is  greater 
than  5.0,  the  explosive-train  safety  will  be  at  least  99-5#. 

18.  The  reliability  of  an  explosive  train  may  exceed  the 
safety  to  such  an  extent  that  the  safety  requirement  Is 
endangered,  or  conversely.  A comparison  of  the  degree  of 
safety,  Dg,  with  the  degree  of  reliability.  Dr,  will  indicate 
whether  the  train  is  over-designed  In  favor  of  either  the 
safety  or  reliability  requirement. 
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50#  Reliable  at  4-1/2  sheets. 

Standard  Deviation  equals  1-1/4  sheets. 


DISTANCE  FROM  ARMED  POSITION  (IN  32NDS  OF  AN  INCH) 

GRAPH  I PROBIT  ANALYSIS  OF  EXPLOSIVE  TRAIN  SAFETY  AND  RELIABILITY 

OF  A GUIDED  MISSILE  FUZE 


DISTANCE  FROM  ARMED  POSITION  (IN  32NDS  OF  AN  INCH) 

GRAPH  3 BRUCETON  ANALYSIS  OF  EXPLOSIVE  TRAIN  SAFETY  AND  RELIABILITY 

OF  A BOMB  FUZE 


DISTANCE  FROM  ARMEO  POSITION  {IN  32NDS  OF  AN  INCH) 

GRAPH  4 BRUCETON  ANALYSIS  OF  EXPLOSIVE  TRAIN  SAFETY  AND  RELIABILITY 

OF  A TORPEDO  FIRING  DEVICE 
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GRAPH  5 BRUCETON  ANALYSIS  OF  EXPLOSIVE  TRAIN  SAFETY  AND  RELIABILITY 

OF  A TORPEDO  ARMINfG  DEVICE 
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GRAPH  6 BRUCETON  ANALYSIS  OF  EXPLOSIVE  TRAIN  SAFETY  AND  RELIABILITY 

OF  A GUIDED  MISSILE  FUZE 


SHEETS  OF  PLASTIC 

GRAPH  7 PROBIT  ANALYSIS  OF  PRIMACORD  PROAGATION 
ACROSS  PLASTIC  SHEET  BARRIERS 
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APPENDIX  A 


THE  BRUCETON  METHOD  OF  SENSITIVITY  TESTING 


INTRODUCTION 

1.  In  order  to  illustrate  the  basic  statistical  application 
of  the  Bruceton  Method,  an  example  of  an  out-of-line  safety 
device  arrangement,  consisting  of  a primer,  a detonator  housed 
in  a slider,  and  a lead,  will  be  used.  Given  a sample  size 
of  30  units  with  a distance  of  one  inch  between  the  "armed" 
and  "safe"  positions,  the  following  procedure  may  be  used  to 
determine  reliability  estimates. 


Test  Procedure 


2.  Choose  a distance,  k,  measured  from  "armed"  position,  at 
which  the  slider  will  be  placed  for  the  first  sample.  This 
distance  should  be  the  position  at  which  roughly  1/2  of  the 
samples  will  fire. 

3.  If  the  first  unit  fires  when  tested  at  a distance  k.  the 
second  one  should  be  tested  at  a distance  k + d,  where  d is  a 
preassigned  increment  between  test  distances.  If  the  first 
unit  does  not  fire,  the  second  one  should  be  tested  at  a dis- 
tance k - d.  Each  succeeding  unit  will  be  tested  at  a 
distance  dependent  on  the  firing  behavior  of  the  previous 
unit:  if  the  previous  unit  fired,  the  next  test  distance 
will  be  d greater  than  the  previous  distance;  or,  if  the 
previous  unit  failed  to  fire,  the  next  test  distance  will  be 
d les3  than  the  previous  distance.  Thus,  the  test  results 
will  be  a sequence  of  fires  (X)  and  failures  (0)  centered 
about  the  50^-firing  distance,  3cR. 

4.  For  the  example  described,  a distance  of  10/32  of  an  inch 
was  chosen  for  the  first  unit  and  Increments  of  1/32  of  an 
inch  for  d.  The  test  results  formed  the  following  sequence: 
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Distance  From 

Armed  Position 

(in  32nd 's  of  an  inch) 

k - 2d  » 8 X 

k - d » 9 *S--X-Q-X-X"X  -X X-X X-X 

k » 10  X-  0 0 0 0 O-O-X  0 0 X 9-9- 

k + d *»  11  — 0 0 X-0 

k + 2d  = 12 0 — 


Calculation  of  the  50#-Rellable  Distance 

5.  The  data  are  tabulated  by  recording  the  number  of  X's  and 
the  number  of  0's  for  each  test  distance.  The  test  distances 
are  assigned  levels,  i (of  equal  unit  increments),  with  the 
level  at  which  no  failures  were  obtained  set  equal  to  zero. 
(Caution;  If  no  failures  were  obtained  at  more  than  one  level, 
the  i = o level  Is  that  distance  which  is  only  d less  than  the7 
level  at  which  failures  were  obtained.  Data  from  test  distances 
less  than  the  i = o level  should  be  disregarded.)  The  calcu- 
lations of  Xp  and  or  Involve  only  the  zeros  or  only  the  X's; 
l.e.,  the  set  which  contains  the  fewer  number  of  results,  n, 
is  used.  The  tabulation  for  the  example  cited  is  as  follows; 


Distance  From 

Armed  Position 

(32nd1 a of  an  inch)  i 

8 0 

9 1 

10  2 

11  3 

12  4 


0 i nA  l2  n1 

0 0 0 

1 10  10 

10  6 12 

3 3 9 

10  0 
n-15  15  A-19  B=31 


»i 

v 

1 

10 

3 

1 

0 


2 


jiA T'.trZ  7.  2 > C 7 ' 1 

unc-sr  I ' = ir_d  7 -racer  cf  J 5 *ec-e  e^ual.  tr.e  i $ 


= 3 + 1.3  = ?.$  (in  32nd1  s of  an  inch) 


Calculation  of  the  Standard  Deviation  (oR) 

7.  The  standard  deviation  is  determined  by  the  following 
8 tops  . 

a.  Find  a value  M * {nB  - Ac)/n^. 

b.  Find  the  value  of  s corresponding  to  the  calculated 
value  M.  Table  1 or  Graph  1 give  s for  values  of  M > .40, 
or  graph  2 gives  s for  values  of  M < .40. 

c.  Determine  aR  = sd . 

For  the  -xfu-ple  the  calculations  are  as  follows: 

M (15  x 31  - (19)2)/(15)2 
« (465  ~ 36l)/225  - 0.46,  then 
from  table  1, 

s = .79,  and 

o_  = .79  (1)  - 0.8  (in  32nd 's  of  an  inch). 
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Calculation  of  Percentage -Reliable  Curve 

8.  Distances  at  which  a given  percentage  of  the  devices  will 
function  may  be  estimated  by  adding  to  or  subtracting  from 
XR  certain  multiples  of  oR.  Percentage-reliable  distances  may 
be  estimated  -'rom  ± toR  by  choosing  the  proper  value  of  t 
shown  below  for  specified  percent  points.  The  percentage 
distances  less  than  50  are  obtained  by  using  the  plus  3ign 
in  the  above  expression,  and  those  greater  than  50  (given  in 
the  parentheses  below)  by  using  the  minus  sign. 


Percent,  p t 


50 

5° 

0 

25 

75 

0.675 

10 

190 

1.282 

5 

(95 

1.645 

1 

99  j 

9) 

2.326 

0.1 

99- 

3.090 

0.01 

[99. 

99) 

3.719 

The  percentage-reliable  curve  for  the  particular  example  cited 
is: 


Percentage - 
Distance  (32nd's 
of  an  inch) 


Percent,  p 

± toR 

± to 

99.99 

-2.98 

6.8 

99.9 

-2.47 

7.3 

99 

-1.86 

7.9 

95 

-1.32 

6.5 

90 

-1.03 

8.8 

75 

-O.54 

9.3 

50 

6 

9.8 

25 

+0.54 

10.3 

10 

+1.03 

10.8 

5 

+1.32 

11.1 

1 

+1 .86 

11.7 

0.1 

+2.47 

12.3 

0.01 

+2  ..98 

12.8 

4 
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Sampling  Error,  ox,  for  the  50/6-Reliable  Estimate 


9»  The  estimate  of  XR  is  subject  to  sampling  error  which  may 
be  oa2 cul- ted  from  oR  by  the  formula, 

ax  = ^/(n)1/2, 

generally  referred  to  as  the  standard  error  of  XR,  where  n is 
the  number  of  0's  (or  X's)  used  in  the  calculation  of  XR  and 
or,  and  G is  read  from  Graph  3 at  the  value  3 determined  in 
computing  or.  For  the  example  the  error  of  XR  is, 

o-  = (0.8)(1.04)/(15)V2 

- 0.22 


Sampling  Error  (aa)  of  the  Standard  Deviation 


10.  The  estimate  of  or  is  also  subject  to  sampling  error, 
which  is  measured  by  its  standard  error,  o0,  given  by  the 
relation 

oa  = ffRH/(n)  ^ * 

where  H is  read  from  Graph  4 using  the  value  of  s determined 
in  computing  oR.  This  measure  calculated  for  the  example  is 

o0  = 0.8(1. 31)/(15)l/2 


Confidence  Intervals  Determined  From  o^  and  aa 


11.  The  standard  errors,  a*  and  a0,  provide  measures  of  pre- 
cision of  the  Xr  estimate  and  the  or  estimate.  Since  XR  and 
or  are  only  estimates  of  their  true  values,  confidence  limits 
for  these  estimates  which  include  the  true  value  will  indicate 
the  precision  of  uhe  estimates.  If  y denotes  an  estimate  and 
Oy  its  standard  error,  confidence  limits  may  be  obtained  by 

y ± toy. 


3 
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where  t Is  a value  given  in  table  2 using  N = (n  - 1)  and 
P = 1 minus  the  level  of  confidence . For  the  illustrated 
example,  the  90#  confidence  limits  of  and  oR  may  be  cal- 
culated as  follows:  1 “ 


N = 15  - 1 » 14, 

P . 1 - .90  « .10, 

which  gives  a t = 1.76; 

therefore,  the  limits  for  Xr  are 

XR  ± 1.76a*  = (9.8)  ± 1.76(0.22) 

= 9*4  to  10.2  (in  32nd' s of  an  inch) 

and  the  limits  for  aR  are 

oR  ± 1.76oa  = (0.8)  ± 1.76  (0.27) 

= 0.3  to  1.3  (in  32nd1 s of  an  inch). 

The  probability  is  10#  that  a value  lies  outside  of  its 
90#  confidence  limits.  Since  this  10#  probability  is  the 
sura  of  the  5#  probability  that  a value  lies  above  its 
greater  90#  confidence  limit  and  the  5#  probability  that 
a value  lies  below  its  smaller  90#  confidence  limit,  it 
can  be  seen  that  the  probability  is  95#  that  a value  will 
lie  above  its  smaller  90#  confidence  limit.  In  other 
words,  a smaller  90#  confidence  limit  is  a 95#  lower 
limit  since  the  probability  is  95#  that  the  true  value 

i i _ _ « u j iL  t ~ w.  ~ — a A/t/  . ...  j - . . - - 

jlaco  auuvc  xu.  oj.miJ.arij , ujic  larger  yoyo  cum xuenoe 

limit  is  a 95#  upper  limit. 


Conclusions  Concerning  the  Degree  of  Reliability 

12.  When  "Qly  ( samples  are  used  to  estimate  the 

reliability  of  . > e-oc . . ive  train,  the  above  confidence 
limits  may  be  useu  to  determine  the  "most  pessimistic"  Dr 
to  be  expected.  By  using  the  lower  limit  for  Xr  and  the’ 
upper  limit  of  oR  to  calculate  D^,  the  "most  pessimistic" 
value  will  be  obtained . The  example  used  gave  a lower  limit 
for  Xr  equal  to  9.4  with  1.3  as  the  upper  limit  for  aR.  The 
"most  pessimistic"  degree  of  reliability  is,  then, 
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9. 4/1. 3 = 7.2, 

as  compared  to  Its  average  degree  of  reliability, 

9. 8/0. 8 = 12.2. 

If  the  "most  pessimistic"  DR  is  greater  than  the  normal 
deviate  of  2.58  which  corresponds  to  99*5$,  then  the  99.5$ 
reliability  is  insured. 


The  Previous  Methods  Applied  to  Safety  Tests 

13.  The  Bruce ton  Method,  as  outlined,  may  be  readily  adapted 
to  safety  tests,  by  recording  "unsafe'1  and  "safe”  results 
rather  than  fires  and  misfires.  If  the  distances  are  measured 
from  "armed"  position  for  the  safety  tests,  the  method  of  cal- 
culation will  be  identical  to  the  preceding  example.  But  to 
obtain  Do,  Xy  must  be  converted  into  units  from  "safe"  posi- 
tion. Also,  the  "most  pessimistic"  D would  be  the  upper 
limit  of  Xy  converted  into  units  from^safe  position,  divided 
by  the  upper  limit  of  oTj. 


7 


TABLE 


NAVJRD  Rgport  2101 


LA  A CO  Ot— rHCMOt— CMVOONrHOJtOCOCOCOCMrHOONC— ACOOOO 
VO  O OJ  AVO  COONOOrHrHrHCMKVJCMOJCMCMCMCMCMrHrHi-.  HriO 
t— 4 O O CM  CO  4 rH  SfOOUAH  t— (OONArHC— (OONArH  t—  l TOMA 

md co  o i-h  co4 >ooo  c^h  t-a>ow4-  m t— oo  oom avo oo 


CM  CM  CVI  CM  (CM  CM  CO  OA  ( 


I CO  (04  4 4 4 4 4 


t-CUVO  : 

ow  mo  : 

lA  CM  CO  : 
MO  CO  ONr 


HO  O rH  O' 
lOOimr-' 

■Ht^-ono' 

I (04  VO  t- 


rH  LACO  O C 
LA  lA  lAMO  « 
CrH  t — CO  ON  L 
\rH  CM4  At 


to co  co  co  cm  «— i o on  t—  aco  o oo 

VO  VO  VO  VO  O VO  VO  A A A A A4 
\HSfOa\AHNrOO\lfiHNtO 
-ONO  CM  CO  A C — 00  O -H  no  AMOOO 


: r-t  i— i r-i  t—t  l-H  CM  CM  CM  CM  (CM  CM  CO  CO  (OKO  CO  <04  4 |4  4 4 4 


VO  ON  CO 
CM  t-O 
M’  O t — 
MO  CO  ON 


>4  ono  on  - 
I 4 IASS 
1 On  A rH  t- 
i CM4vo  c— < 


rH  ACO  ok 
ONONONOK 
ION  AH0O|= 
v O CM  4 At 


co  co  ro  ro : 
o o o o c 
OVO  CM 00 
- ON  o CM  CO 


rHOONt— ACOrHOO 
O O ON  ON  ON  ON  ON  00 
OVOHNCOAAH 
\ C— CO  O rH  CO-4-  VO  00 


rH  rH  pH  rH  |rH  CM  CM  CM  CM  (CM  CM  COCO  COlCO  CO  (04  4 (4  4 4 4 


CO  C—  rH  : 
ArH-4-  . 
CM  ON  Ar 
VO  f— ON> 


■ CO  t — ONCO  A O A 00  O 
) CO  ON  O rH  CM  CO  CO  CO  4 
lf-(OOV0CMC04OV0 
I CM  -4-  MD  C — ON  O CM  -4"  A 


I CO  CO  CO  COt 
• .4-  -4-  4 4 = 
iOO-4-  OMO 
-00  O CM  CO 


1 rH  O ON  C — A CO  rH  ON 
-44CO(OCOCOCOCM 
i 00  4 O VO  CM  CO  -4"  O 
NVOCO  O rH  CO-4"  VO  CO 


CM  CM  CM  CM  )0J  CM  COCO COk 


14  4- 14-  -4-  -4"  4 


t- — 4"  ONr 
r-AM 
O £—  tok 
VO  t—  ONr 


I rHVOOO  t—  = 
• CM  CO-4"  A . 
HO  CM  CO  -4"  i 
I CM  4 At-  t 


• 0400  O 

> t—  t— t— ao 

) VO  CM  004 
NO  CM  CO  A 


I co  co  co  cok 

100  00  oo  oot 

l VO  CM  CO  4 S 
-CO  O rH  <04 


1 CM  O ON  t— lA  CO  rH  ON 

> co  oo  t-  t— b—  t—  t— vo 
! 'O  OJ  00  -4-  p VO  OJ  CO 
NVOCO  0NrnkO4V0  t— 


iHrHrHrHkHCMCMCMCMlCMCM  COCO  COlCO  CO  CO  (04  |4 4 4 4 


C rH  NO  C 
ON  ON  rH  ( 
00  ACM 
At- ON 


N ON  A E — A 
1 A t-CO  ON 
>4  OVO  CM 
I CM  4 At— 


ON4  t— Ok 
OHHCMk 

1 AH  t-cok 
I O CM  CO  aL 


i CO  CO  COCO 
I CM  CM  CM  CM 
NAH  t— CO 
ICO  o rH  CO 


NCM  rH  OOO 
I CM  CM  CM  rH 
NAH  t—  CO 
VO  CO  ON  rH 


14  rH  ON 
I rH  rH  O 
NAH  t— 
I 4 VO  t— 


rH  rH  (Hi 


I CM  CM  CM  CM  I 


) ro  coko  CO  CO  (04  |4 4 4 4 


rH  t— 4 
rH  CM  A 
t— 4 O 
AC— ON 


-C—  COVO  AVO  ON  CO  t—  On! 


- ON  rH  CM  COf4  4 A A At 
'CM  ON  A rH  F—  CO  ON  A rH  I 


VM  M/\U  li — if.'-  CO  ON  A rH  | 
CMCOAC— nOrHCOAl 

rCHHHHCMCMCMCM 


I CM  (OCOfOC 

I VO  'O  VO  VO  - 
-COON  AH  t 
ICO  ON  rH  CO: 


NCM  rH  OOO  - 
I VO  VO  VO  Al 
-COON  AH  [ 
VO  C — ON  rH  ( 


14  CM  ON 
N A A4 
-COON  A 
14  At— 


CM  CM  CO  (OVO  CO  CO  (04  4 4 4 4 


CO  (OCM  AACM  AA( 
rH  VO  ON  rH  (O  AVO  t—  : 
ACM  CO  AH  C—  CO  ON  l 
A t— 00  O CM  CO  AVO  : 


ICO  COC— ON 
ICO  ON  ON  ON 
N rH  t— COON 
I O rH  (04 


I CM  CO  CO  CO 
! O O O O 
I CM  004  O 
ICO  ON  rH  CO 


NCM  rH  OCO 
I O O O ON 
I CM 004  ON 
•VO  t— ONO 


>4  CM  O 
NON  ON  ON 
N rH  C— CO 
14  AC— 


r-\  r-l  r-\  r-i\ 


icuajcuojlojcvajpo  cok 


O CO  (04  J4  4 4 4 


VO  00  0NKO4  04  4 
rH  ON  CM  At— ONO  rH 
COO  t— 1(0  ON  A CM  00 
A C— 00  p H CO  AVO 


I CO  COVO  ON  r 
I CM  CO  CO  CO  : 
• O VO  CM  00  : 
I O rH  (04  . 


1 CM  rococo 
■4444 
OVO  CM  CO 
I 00  ON  rH  CM 


NCM  rH  OOO 
-444  CO 
- OVO  CM  CO 
-VO  t— OnO 


54  CM  O 
N CO  (OCO 
r OVO  CM 
J4  At— 


CM  CM  CM  CM  KM  CM  CM  (O  (OKO  CO  CO  (04  {4444 


CM  CM  t— P CM  ON  CO  CO 
O COVO  (ON  rH  OJ  4 A 
rH  ON  A|H  CQ  4 O VO 
AVO  00  p-r  H fO  uOVO 


I fj  ri(*  I U I' 

> I • • • « 

jr-t  rH  rH  rH  1 


I t—  CM  VO  ONr 
) VO  t—  C-  t-t 
I CO  4 ovo  k 
I ON  H 004  t- 


1 CM  ro  co  rok 

I CO  CO  00  cot 

ioo4  ovok 

1 C—  ON  rH  CM  r 


NCM  rH  OOO  t — 4 CM  OCO 
1 00  GO  GO  t— t — t— t—  t— VO 
1 00  4 O VO  CM  00  4 O VO 
■ A t—  ON  O “CM  CO  A t-CO 


CM  CM  CM 


CM  CM  CO  (OkO  CO  CO  (04  444  44 


OOO  jo  OOOO  k-H  CM  (04  A O C— CO  ON  O rH  CM  (04  A O t— 00  ON  O 
(04  A|JD  t-CO  ON^[_^  r-4  H rH  hH  rH  rH  rH  CU*|cM*  CM*  CM*  CM*  CM  (CM*  CM*  CM*  CM*  CO 


I 


NAVORD  Report  2101 


TABLE  2 

Values  of  t for  N = (n  - 1) 
versus  P = 1 - level  of  confidence 

(P) 


N 

.90 

.70 

♦ 50 

.30 

.10 

.05 

.02 

.01 

1 

.16 

.51 

1.00 

1.96 

6.31 

12.7 

31.8 

63.7 

2 

.14 

„44 

.82 

1.39 

2.92 

4.30 

6.96 

9.92 

3 

.14 

.42 

.76 

1.25 

2.35 

3.18 

4.54 

5.84 

4 

.13 

.41 

. 7b 

1.19 

2.13 

2.78 

3^75 

4.60 

5 

.13 

.41 

.73 

1.16 

2.02 

2.57 

3.36 

4.03 

6 

.13 

.40 

• 72 

1.13 

1.94 

2.43 

3.14 

3.71 

7 

.13 

.40 

.71 

1.12 

1.90 

2.36 

3.00 

3.50 

8 

.13 

.40 

.71 

1.11 

1.86 

2.31 

2.90 

3.36 

9 

.13 

.40 

.70 

1.10 

1.83 

2.26 

2.82 

3.25 

10 

.13 

.40 

.70 

1.09 

1.8l 

2.23 

2.76 

3.17 

12 

.13 

.40 

.70 

1.08 

1.78 

2.18 

2.68 

3.06 

14 

.13 

-39 

.69 

1.08 

1.76 

2.14 

2.62 

2.98 

16 

.13 

.39 

.69 

1.07 

1.75 

2.12 

2.58 

2.92 

18 

.13 

.39 

.69 

1.07 

1.73 

2.10 

2.55 

2.88 

20 

.13 

.39 

.69 

1.06 

1.72 

2.09 

2.53 

2.84 

22 

.13 

.39 

.69 

1.06 

1.72 

2.07 

2.51 

2.82 

24 

.13 

.39 

.68 

1.06 

1.71 

2.06 

2.49 

2.80 

26 

.•13 

.39 

.68 

1.06 

1.71 

2.06 

2.48 

2.78 

28 

.13 

.39 

.68 

1.06 

1.70 

2.05 

2.47 

2.76 

30 

.13 

.39 

.68 

1.06 

1.70 

2.04- 

2.46 

2.75 

40 

.13 

.39 

.68 

1.05 

1.68 

2.02 

2.42 

2.70 

60 

.13 

.39 

.68 

1.05 

1.67 

2.00 

2.39 

2.66 

120 

.13 

.39 

.68 

1.04 

1.66 

1.98 

2.36 

2.62 

00 

.126 

.385 

.674 

1.036 

1.643 

1.960 

2.326 

2.576 

GRAPH  I CURVES  FOR  ESTIMATING  THE  STANDARD  DEVIATION  (M>0.40) 


GRAPH  2 CURVES  FOR  ESTIMATING  THE  STANDARD  DEVIATION  (M<0.40) 
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GRAPH  3 CURVE  FOR  ESTIMATING  THE  STANDARD  ERROR  OF  THE  MEAN  (s  > 0.5) 
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GRAPH  4 CURVE  FOR  ESTIMATING  THE  STANDARD  ERROR  OF  a (s>0.5) 
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GRAPH  5 CURVE  FOR  ESTIMATING  STANDARD  ERRORS  OF  THE  ESTIMATES 

OF  THE  MEAN  AND  STANDARD  DEVIATION 
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APPENDIX  B 


THE  PROBIT  METHOD  OF  SENSITIVITY  TESTING 


INTRODUCTION 

1.  In  order  to  illustrate  the  Probit  Method,  consider  a 
safety  test  on  an  explosive  train  which  contains  a primer,  a 
slider  detonator,  and  a lead.  The  slider  safety  device  travels 
a distance  of  16/32  of  an  inch  between  "armed"  and  "safe" 
positions.  Assume  that,  when  the  slider  detonator  is  in  any 
position  in  which  it  can  be  initiated  by  the  primer,  it  will 
always  initiate  the  lead.  Thus,  the  safety  of  the  train 
is  dependent  on  the  slider-detonator  to  lead  transfer,  since 
the  most  unsafe  condition  exists  at  this  transfer.  Therefore, 
estimates  of  safety  will  be  obtained  by  purposely  initiating 
the  slider  detonator  at  greater  distances  than  the  primer 
transfer  allows  and  recording  the  effect  on  the  lead  as  "safe" 
(0)  or  "unsafe"  (X) . In  this  manner,  if  spurious  initiation 
of  the  slider  detonator  occurs,  it  may  be  determined  whether 
the  safety  of  the  train  is  adequate  or  not,  as  the  case  may  be. 


Test  Procedure  for  the  Probit  Method 


2.  Determine,  roughly,  the  distance  from  armed  position 
which  produces  an  unsafe  condition  about  50%  of  the  time 
when  the  slider  detonator  is  purposely  Initiated.  The  50% 
distance  may  be  roughly  estimated  by  testing  a small  sample 
of  5 by  means  of  the  up-and-down  procedure.  For  illustra- 
icn,  the  following  results  were  obtained  from  an  actual  tests 


Distance  from 
Armed  Position 


* 


■e X 


5/32 

6/32 

7/32 


1 
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3.  Samples  of  at  least  10  each,,  including  the  first  3 tested 
above,  should  be  tested  at  distances  around  the  50$  distance, 
thus  obtaining  a range  of  percentage -distances . Such  percent- 
age distance.?,  may  be  obtained  in  any  order  and  do  not  require 
an  equal  niafvcr  of  samples  for  each  distance  tested;  nor  do 
the  distances  tested  need  to  be  equally  spaced.  The  high  and 
low  percentage  distances  should  be  estimated  by  larger  samples 
than  those  near  the  50$  distance,  since  they  are  less  accurate 
when  estimated  from  small  samples.  The  results  obtained  for 
the  illustrative  example  are  tabulated  as  follows: 

Distance  from 

Armed  Position  Percent 


(in  32nd* 3 of  an 
inch) 

Unsafe 

X 

Safe 

0 

Unsafe 

P 

5 

XXXXOXXXXXXXXXXXXX 

17 

1 

94.7 

5-3/4 

XXXXX0X0X0 

7 

3 

70 

6 

0XXXXX0XXX 

8 

2 

80 

6-1/4 

oxoooxooox 

3 

7 

30 

7 

000000X000 

1 

9 

10 

Total 

w~ 

22" 

Description  of  a 

"Probit" 

4.  The  Probit  Analysis  Is  based  on  the  transformation  of  the 
percentage  unsafe  into  "Probits".  To  illustrate  this  step 
clearly,  some  understanding  of  the  theoretical  normal  curve 
is  necessary.  The  theoretical  normal  curve  i£  determined  by 
its  two  parameters,  the  50$  unsafe  distance,  Xy,  and  the 
standard  deviation,  oy.  Any  specified  percentage -unsafe 
distance  may  be  obtained  by  integration  of  the  normal  curve 
equation.  Tables  of  this  Integration  are  available  in  most 
statistical  texts  which  give  the  percentage  of  the  area  under 
the  curve  for  various  distances,  x,  measured  from  in  terms 
of  ouj  00  defined  as, 

. _ * “ Xu 

au 
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Thus,  any  desired  percentage -unsafe  distance,  x,  may  be  __ 
obtained  from  Xy  and  ay  by  x = Xy  z tay . Now  if  ay  and  Xy 
are  not  known  but,  say,  a 30%  distance  is  obtained,  the  table 
of  the  nc>".c-a1  curve  areas  will  give  a t value  that  represent.3 
the  difference  between  Xy  and  the  30%  distance  in  terms  of 
dy.  Thus,  the  "Probit"  of  a percentage-unsafe  distance,  x, 
is  the  number  of  standard  deviations  that  x is  removed  from 
the  50%  distance,  and,  for  simplicity  in  calculation,  the 
probit  of  Xy  is  designated  as  5.0.  For  percentages  below 
50%  the  probits  are  determined  by  subtracting  the  number  of 
standard  deviations  from  5.0,  and  for  percentages  above  50%, 
the  number  of  standard  deviations  are  added  to  5-0.  A probit, 
Y,  may  be  defined  as, 

Y = 5 + ~ *u) . 

°U 

Table  1 presents  a complete  tabulation  of  the  transformation 
of  percentages  to  probits. 


Analysis  of  Test  Data  by  the  Provisional  Problt.  Line 


5.  Step  1 of  the  analysis  is  the  transformation  of  the 
observed  percentages  to  probits,  designated  as  "empirical 
probits".  By  using  Table  1,  these  values  may  be  read  directly. 

6.  Step  2 involves  plotting  a graph  of  the  empirical  problts 
as  a function  of  the  distances  tested.  A provisional  straight 
line  is  then  drawn  that  best  fits  the  empirical  probits, 
placing  the  line  by  eye.  For  each  of  the  distances  tested, 
the  value  of  the  ordinate  of  the  provisional  line  is  read. 

These  are  termed  tne  "provisional  problts",  Y. 

7 . If  the  provisional  line  is  drawn  so  the  difference  between 
the  empirical  probits  and  the  provisional  problts  is  a minimum, 
and  no  other  line  will  fit  the  empirical  problts  more  accu- 
rately, the  more  precise  mathematical  procedures  for  deter- 
mining the  problt  line  may  not  be  needed. 

8.  The  equation  for  the  provisional  line  may  be  obtained 
directly  from  the  graph.  The  estimate  of  Xy  is  the  x value 
which  corresponds  to  the  problt  5.0.  An  estimate  of  ay  is 
the  increase  in  x for  a unit  increase  of  the  problt  values. 
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The  slope  of  the  line,  b,  is  given  by  l/oy.  The  equation  of 
the  line  can  then  be  obtained  by  substituting  the  values  of 
X^j  and  b in  the  equation 


Y = 5 + b(x  - Xu) . 

9.  As  an  illustration,  the  test  results  presented  in  para- 
graph 3 may  be  analyzed  by  the  methods  of  paragraphs  5 
through  8 as  follows: 

Distance  from 
Armed  Posicion 


(in  32nd' s of 
an  inch) 

X 

Number 

Tested 

n 

Percentage 

Unsafe 

P 

Empirical 

Problts 

Provisional 

Probits 

Y 

5 

18 

94.7 

6.62 

6.6 

5-3/4 

10 

70 

5.52 

5.5 

6 

10 

80 

5.84 

5.2 

6-1/4 

10 

30 

4.48 

4.8 

7 

10 

10 

3.72 

3-7 

The  empirical  problts  were  plotted  in  Graph  1 and  the  provi- 
sional line  estimated.  The  provisional  probits  were  obtained 
from  the  provisional  line  of  Graph  The  Xy  is  the  x value 
corresponding  to  the  probit  5. 0 , which  is  6.12  measured  in 
32nd’s  of  an  inch.  A unit  Increase  in  the  probit  value  of 
the  provisional  line  provides  a decrease  in  the  value  of  x 
equal  to  0.69,  measured  in  32nd1 s of  an  inch.  This  value  is 
an  estimate  of  du>  hut  since  it .represents  a decrease  in  x 
rather  than  an  increase,  it  must  be  considered  as  a negative 
value  in  determining  b.  Thus,  the  slope,  b,  of  the  provi- 
sional line  i£  l/ cry  or  1/-Q.69,  which  is  equal  to  -1*45.  By 
substituting  Xu  = 6.12,  and  b = -1.45  in  the  equation 

Y « 5 + b(x  - Xu), 

one  may  obtain  the  equation  of  the  provisional  probit  line 
in  terms  of  x,  or  distance  from  armed  position.  Thus, 
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Y = 5 - 1.45(x  - 6.12) 
or,  Y = 13.87  - 1.45.x 

10.  The  equation  of  the  provisional  line  provides  a means 
for  estimating  the  distance  at  which  a specified  expected 
probit  occurs.  If  the  distances  related  to  problts  between 
3.0  and  7.0  are  calculated,  the  percentages  corresponding 
to  these  problts  may  be  obtained  from  Table  1,  and  then  the 
relationship  of  distances  versus  percentage-unsaf e may  be 
plotted . 

11.  The  foregoing  method  does  not  provide  a mear^s  for  deter- 
mining the  error  of  the  Xu  or  ajj,  so  it  should  not  be  used  as 
the  complete  analysis  unless  D3  is  considerably  larger  than 
5.0.  Dg  is  obtained  by  dividing  the  Xu,  measured  from  safe 
position,  by  au  obtained  from  the  provisional  line.  £or  the 
illustrative  example: 

safe_posltion  is  16/32  of  an  inch  from  armed  position, 
and  Xy  is  6.12/32  of  an  inch  from  armed  position;  there- 
fore,  t r r 10 

DS  - 0.69  - 9- 88/. 69  = 14.3 

wiiich  indicates  that  the  safety  is  more  than  adequate  and 
that  further  calculation  is  not  necessary.  But  for  illus- 
tration and  comparison,  this  example  will  be  used  to  demon- 
strate the  calculation  of  the  Probit  Line. 


The  Method  of  Calculating  the  Probit  Line 

12.  The  method  of  calculating  the  Problt  Line,  based  on  the 
maximum  likelihood  theory  of  reference  (c),  has  been  simpli- 
fied considerably  by  use  of  tables  which  are  appended  to  this 
report.  The  explanation  of  the  theoretical  basis  of  the 
method  has  been  omitted,  so  the  following  example  is  only  an 
illustration  of  the  calculations  and  application  of  the  tables. 

13.  The  preceding  tabulation  in  paragraph  9 for  estimating 
the  provisional  line  Is  a necessary  part  of  che  calculations 
for  the  probit  line.  Thus,  steps  1,  2,  and  3 are:  (1)  trans- 
forming the  percentages  to  empirical  probits  by  means  of 
Table  1,  (2)  plotting  the  empirical  probits  and  estimating 


5 


NAVORD  Report  2101 


the  provisional  probit  line  by  eye,  and,  (3)  reading  the 
provisional  probits,  Y,  from  the  provisional  line.  Step  4 
incorporates  the  weighting  coefficient,  w,  which  is  read 
from  Table  2 for  each  provisional  probit  value.  The  total 
number  of  trials  at  each  distance  is  n,  which  is  multiplied 
by  w and  entered  in  the  tabulation  under  nw.  Step  5 deter- 
mines a "working  probit",  y,  for  each  empirical  percentage 
versus  its  provisional  probit.  This  value  is  read  from 
Table  3.  For  the  example  cited,  the  tabulation  of  the  data 
is  as  follows: 

Distance  from 

Armed  Position  Number  Percent  Empirical  Provi-  Weighting  Working 
(32nd 's  of  an  Tested  Unsafe  Probits  sional  Factor  Probits 

inch)  Probits 


X 

n 

P .... 

Y 

nw 

y 

5 

- r\ 

10 

94.7 

6.62 

6.6 

4.28 

6.62 

5-75 

10 

70 

5.52 

5.5 

5.81 

5.52 

6 

10 

80 

5.84 

5.2 

6.27 

5.76 

6.25 

10 

30 

4.48 

4.8 

6.27 

4.49 

7 

10 

10 

3.72 

3-7 

Znw 

3.36 

= 

3.72 

14.  The  products  nwx  and  nwy  are  tabulated  and  added;  division 
of  the  totals  by  Znw  gives  the  means,  x and  y.  The  values  nwx 
are  multiplied  by  x and  summed  to  give  Znwx  . Also,  they  are 
multiplied  by  y and  added  to  give  Znwxy.  The  values  nwy  are 
multiplied  by  y and  added  to  give  Znwy2.  The  sums  of  squares 
and  products  about  the  mean  are  then  obtained  as  follows: 


(a) 

Sxx  = Znwx2  - 

( Znwx)2/Znw, 

(b) 

Sxy  = Znwxy  - 

( Znwx ) { Zr  wy ) /Znw 

(c) 

Syy  = Znwy2  - 

(Znwy)2/Znw. 

15.  The  equation  of  the  Probit.  Line  may_be  obtained  by  sub- 
stituting in  the  equation  Y = y + b(x  - x)  the  calculated 
values 
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y 


Znwy 
Znw  ’ 


■7  _ Znwx 

~ 2H^T' 

and  b = Sxy/SXx>  which  give  a linear  rela- 
fc  1 onn’'-'.p  between  Y,  in  problts,  and  x,  distance  from  armed 
posxtlon. 


16.  The  previous  example  provides  the  following  calculations 
as  Illustration  of  paragraphs  14  and  15: 


nwx 

21.4000 

33-4075 

37-6200 

39.1875 

23.5200 


nwy 

23.3336 

32.0712 

36.1152 

28.1523 

12.4992 


nwx 


107.000000 

192.093125 

225.720000 

244.921875 

164.640000 


nwy*- 

187.568432 

177.033024 

208.023552 

126.403827 

46.497024 


nwxy 

141.668000 

184.409400 

216.691200 

175.951875 

87.494400 


=137  :i'7l5Ihwx2=9'3'4  .'3'7500'0Znwy2=?4'575"?55?9Znwxy=ro8';  214875 


Znwx2  = 934.375000 
Znwy2  = 745.525859 
Znwxy  = 806.214875 


( Znwx) (Znwy)/Znw  = 818.78032 
(Znwx)2/Znw  = 926.00493 
(Znwy)2/Znw  = 723-97154 


thus, 

(a)  Sxx  = 934.375000  - 926.00493  = 8.37007, 

(b)  Syy  = 745.525859  - 723.97154  = 21.55432,  and 

(c)  Sxy  = 806.214875  - 818.78032  = 12.56544. 

The  preceding  measures  allow  the  calculations  of  b,  x and  y 
for  substitution  in  the  equation  for  the  Probit  Line, 

Y = y + b(x  - x) . The  calculated  values  for  the  example  are 


and. 

Therefore, 


7 = 137.1715/25.99  = 5.278 
x = I55.I350/25.99  = 5.969 
b = -12.56544/8.37007  = -1.5012. 

the  Probit  Line  equation  for  the  example  is 
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Y = 5.278  - 1.501(x  - 5.969) 
or  Y = 14.237  - 1.501x 

Calculation  or  the  50^-Unsafe  Distance,  X]j,  and  Standard 
Deviation,  ctu 

17.  The  estimated  Xy  is  obtained  by  substituting  in  the 
Problt  Line  equation  the  value  5.0  for  the  probit  value  Y and 
solving  for  x.  For  the  example 

5.0  = 14.237  - 1.501x 

or  x = 9.237/1.501  = 6.15. 

Thus  Xy  equals  6.15#  measured  in  32nd1 s of  an  inch.  The 
estimated  au  is  equal  to  1/b,  or  for  the  example, 

1/1.501  which  equals  0.67  (in  32na's  of  an  inch). 


The  Degree  of  Safety,  Ds,  and  Comparison  Between  the  Provi- 
sional  Estimates  and  Calculated  Estimates 


18.  The  degree  of  safety,  Dg#  is  obtained  by  dividing  Xy, 
measured  from  safe  position,  by  ou#  obtained  from  the  calcu- 
lated problt  line.  From  paragraph  17,  Xy  for  the  example  is 
6.15,  measured  from  arme.d  position  in  32nd  's  of  an  inch, 
which  gives  9.85  from  safe  position  (16  - 6.15),  and  cy  is 
0.67.  Therefore,  .Dg  is  found  to  be  9.85/O.67  = 14.7. 

19.  The  method  of  calculating  the  problt  line  is  an  appro- 
ximation that  results  in  convergence;  therefore,  if  the 
provisional  line  differs  considerably  from  the  calculated  one, 
it  is  advisable  to  perform  another  cycle  of  computations  and 
compare  the  first  cycle  with  the  second  to  see  if  convergence 
has  been  approximated.  Thus,  a comparison  between  the  eye 
estimates  and  the  calculated  estimates  seems  in  order.  For 
the  example,  the  following  measures  may  be  compared: 
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Provisional  Estimates 

Calculated 

Estimates 

Xv  = 6.12 

Xu  = 

6.15 

cr:  « 0.69 

11 

O 

0.6? 

dt  Line, 

Probit  Line, 

Y = 13.87  - 1.45x 

Y = 

14.24  - 1.50x 

»S  - “-3 

DS  = 

14.7 

The  difference  between  the  two  estimates  is  so  small  that  no 
further  calculation  of  the  Probit  Line  is  needed."  Should 
there  be  any  doubt  concerning  the  convergence  of  the  provi- 
sional line  and  the  calculated  line,  the  calculations  described 
in  paragraphs  20  and  21  may  be  performed  to  determine  whether 
another  cycle  of  computations  is  necessary. 


Sampling  Errors  of  the  50%-Unsafe  Distance  and  Slope,  b,  of 
the  Calculated  Probit  Line 


20.  The  error,  a to  expected  in  the  calculated  estimate 
of  may  be  obtained  by 

1 

°x  = b ( Znw ) * 


or  for  the  example 


= 1. 50(25. 99)1/^ 


1 

1.50(5.10) 


0.13. 


If  the  limits  obtained  from  X..  ± include  the  provisional 
Xu>  a further  cycle  of  computations  is  not  necessary.  But  if 
the  provisional  Xy  is  outside  these  limits  another  cycle  should 
be  performed  and  comparison  made  between  cycle  1 and  cycle  2 
estimates.  The  provisional  X..  of  6.12  is  well  within  the  limits 
of  6.02  and  6.28,  indicating  that  sufficient  convergence  of  the 
% es  timates  has  been  obtained. 

21.  The  error,  o^,  to  be  expected  in  the  calculated  slope,  b. 

Is  determined  from 
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°b  = <“s — ) • 

0 xx 

The  error  for  the  example  is: 

ab  = (H7 h)1/2  = (O.lloA)1/2  = 0,35 

If  the  -limits  obtained  from  b ± ab  include  the  provisional 
estimate  of  b,  a further  cycle  of  computations  Is  not  necessary. 
But  if  these  limits  do  not  include  the  provisional  b,  another 
cycle  should  be  performed.  The  provisional  b of  1.45  for  the 
example  is  well  within  the  limits  of  1.15  and  1.85,  which  indi- 
cates that  sufficient  convergence  has  been  obtained. 

22.  The  90 % confidence  limits  of  Xy  may  be  determined  by 

*U  ± tax> 

where  t is  obtained  from  Table  2 of  Appendix  A and  is  the 
normal  deviate  for  P = .10,  or  t = I.65.  The  example  gives 

6.15  ± 1.65(0.13) 

or,  6.15  ± 0.21,  which  gives  limits  of 
5.94  to  6.36. 

23.  Now,  the  90$  confidence  limits  of  ay,  using  t = I.65, 
are  obtained  by 


or,  for  the  exampl* 


b + tah, 


1.50  ± 1.65(0.35)  = 1.50  ± 0.58, 


which  gives  limits  of  0,92  to  2.08.  By  obtaining  the  recip- 
rocals of  the  confidence  limits  of  b,  one  may  determine  the 
limits  of  <jy,  which  are 

1/0.92  to  1/2.08, 

or,  1,09  to  0.48 

24.  The  "most  £essimlstic"  Ds  may  be  calculated  from  the  90# 
upper  limit  of  Xu,  measured  from  safe  position,  and  the  90^ 
upper  limit  of  cru  as  follows,  where  X^  represents  the  total 
distance  from  arrrri  to  safe  position: 
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_ XT  - (Xu+tax) 

S l/(b  + tab) 

p',r*  the  example,  Xjj  + ta~  Is  6.36,  and  l/(b  + ta^)  Is  1.09; 
therefore,  Dg  is  calculated  as  follows: 

Do  = ---■  6-?6  = 8.8 

s 1.09 


This  Dg  is  considerably  better  than  2.58,  the  normal  deviate 
corresponding  to  99-5$>  therefore,  the  device  would  be  con- 
sidered safe  . 


A Test  for  Goodness  of  Fit 


25.  In  some  cases,  it  may  be  doubtful  whether  the  probit  line 
is  actually  a measure  of  the  relationship  between  percentage- 
unsafe  and  distance  tested.  If  the  empirical  probits  are 
scattered  in  such  a fashion  that  the  probit  line  could  be 
drawn  in  several  positions  by  eye  estimate,  it  is  advisable 
to  calculate  the  probit  line  and  obtain  a measure  of  "good- 
ness of  fit"  termed  the  chi-squared  test. 


26.  A chi-squared  value, ✓X2,  is  calculated  as  follows: 

- 2)  “ Syy  - <Sxy/Sx*) 

where  (k  - 2)  refers  to  the  "number  of  degrees  of  freedom", 
or  the  number  of  test  distances  minus  2.  By  referring  to  a 
chi-squared  table  (available  in  reference  (c)  or  other 
standard  statistical  references),  one  may  ascertain  the  pro- 
bability of  obtaining  such  a value  with  (k  - 2)  degrees  of 
freedom,  providing  the  probit  line  is  the  true  relationship 
between  percentage  unsafe  and  distance  tested.  For  the 
example: 

/xfk  - 2)  = 21 *55^32  - (-12.56544)2/8.37 
= 21.55432  - 18.86383  = 2.69049 


with  (k  - 2)  degrees  of  freedom  equal  to  3«  The  chi-squared 
table  gives  a probability  of  0.47,  which  means  that  practi- 
cally half  the  time  one  would  expect  a larger/X^  value  than 
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2.69  from  samples  explained  by  the  calculated  probit  line. 
In  other  words,  the  px-oblt  line  is  a "good  fit". 

27.  A probit  line  which  gave  a chi-squared  value  corres- 
ponding to  a probability  less  than  .05  is  considered  a poor 
fit  and  would  not  be  acceptable  as  the  true  relationship 
between  percentage  unsafe  and  distance  tested.  For  such 
cases,  one  might  determine  some  function  of  x for  which  a 
calculated  probit  line  would  fit  more  closely,  such  as  log 
x,  x , or  1/x,  etc.  The  log  x transformation  is  used  to 
some  extent;  in  which  case  the  probit  equation  takes  the 
form 


Y = y + b(log  x - log  x), 
and  similarly  for  other  transformations. 
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TABLE  3 

WORKING  PROBITS  (y) 


provisional  probit, ,Y 
2.1  2.2  ! 2.3  I 2.4  I 2.5  2.6  I 2.7 
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TABLE  3 (CONTINUED) 
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% 

3-0 

provisional  probits 

Y 

p 

J.J.4 

! Q 0 

° — 

1 3.3 

3.4 

3-5 

3.6 

3.7 

. 3.8 

3.9  1 

0 

i 

2.579 

i 

2.662 

2.745 

2.826 

2.906 

2.984 

3.061 

3.135 

| 3.207 

i 

3.277 

1 

2.764 

I2.815 

2.872 

2.932 

2.996 

3.061 

3.127 

3.193 

3.259 

3.323 

2 

.949 

.967 

.998 

3.039 

3.086 

.139 

.194 

.252 

.310 

.369 

3 

3.134 

3.120 

3.125 

.145 

.176 

.216 

.261 

.310 

.362 

.415 

4 

• 319 

,272 

.252 

.251 

.267 

.293 

.328 

.36° 

.413 

.461 

5 

.505 

.424 

.378 

.358 

,357 

.37° 

.395 

.427 

,465 

.507 

6 

3.690 

3.577 

3.505 

3.464 

3.447 

3.447 

3.461 

3.485 

3.516 

3.553 

7 

.875 

.723 

.632 

.570 

.537 

.525 

.528 

.544 

.568 

.599 

8 

4,060 

.882 

.758 

.677 

.627 

.602 

• 595 

.602 

.619 

.645 

9 

.246 

4.034 

.885 

.783 

.717 

.679 

.662 

.660 

.671 

.690 

10 

431 

,186 

,4,012. 

.889 

.808 

.756 

.728 

.719 

.722 

.736 

11 

4.616 

4.339 

4.138 

3.996 

3.898 

3.334 

3-795 

3 ' 777 

3-774 

3.782 

12 

.801 

.491 

.265 

4.102 

.988 

.911 

.862 

• 835 

.825 

.828 

13 

.986 

.644 

• 391 

.208 

4.078 

.986 

.929 

894 

.877 

.874 

14 

5.172 

.796 

.518 

.315 

.168 

4.065 

.996' 

.952 

.928 

.920 

15 

.357 

.948 

.645 

.421 

.258 

.142 

4.062 

4.010 

.980 

.966 

lb 

5.542 

5.101 

4.771 

4.527 

4.348 

4.220 

4.129 

4.069 

4.031 

4.012 

17 

.727 

.253 

.898 

.634 

.439 

.297 

.196 

.127 

.083 

.058 

18 

.913 

.406 

5.025 

.740 

.529 

.374 

.263 

.185 

.134 

.104 

19 

6.098 

.558 

.151 

.846 

.619 

.451 

.330 

.244 

.186 

.149 

| 20 

.283 

.710 

.278 

.953 

,709 

.528 

.396 

.302 

.237 

.195 

\2 F 

6.468 

5 -863 

5.405 

5.059 

4.799 

4.606 

4.463 

4.361 

4.289 

4.241 

22 

.653 

6.015 

••531 

.165 

.889 

.683 

• 530 

.419 

.340 

.287 

23 

.839 

.168 

.658 

.272 

..979 

.760 

.597 

.477 

.392 

.333 

24 

7.024 

.320 

.785 

.378 

5.070 

.837 

.664 

.536 

.443 

.379 

25 

.20? 

.472 

.484 

.160 

.914 

.730 

.594 

.495 

.425 

26 

7.394 

b!625: 

s!  038 

5.591 

5.250 

4.992 

4.797 

4.652 

4.p46 

4.471 

27 

.580 

.777 

.165 

.697 

.340 

5.069 

.864 

.711 

• 598 

.517 

28 

.765 

.930 

.291 

.803 

,430 

.146 

.931 

.769 

.649 

.563 

29 

.950 

7.082 

.418 

.910 

.520 

.223 

.997 

.827 

.701 

. 608 

30 

8.135 

.234 

.545 

6.016 

.610 

.300 

5.064 

.886 

.752 

.654 

31 

8.320 

7.387 

6!  671 

6.122' 

5.701 

IF37F 

5.131 

4 1944 

4 .804 

4.700 

32 

.506 

.539 

.798 

.229 

.791 

.455 

.198 

5.002 

.855 

.746 

33 

.691 

.692 

.925 

.335 

.881 

.532 

.265 

.061 

.907 

.792 

34 

.876 

.844 

7.051 

.441 

.971 

.609 

.331 

.119 

.958 

.838 

35 

9.061 

.9?6 

,178 

.548 

6.061 

.687 

.398 

,177 

5.010 

.884 

36 

9^W 

8.149 

7.305 

6.654 

6.151 

5.764 

5.465 

5.236 

5.061 

4.930 

37 

.432 

.301 

.431 

.760 

,242 

.841 

.532 

.294 

1V3 

O 7A 
•/:  w 

38 

.617 

.%4 

.558 

.867 

.332 

.918 

.599 

.353 

.164 

5.022 

39 

.802 

.606 

.685 

.973 

.422 

• 995 

.665 

,4ll 

,216 

.068 

40 

.987 

.758 

.811 

7,07? 

.512 

6.073 

,732 

.469 

,267 

,1.134 

41 

8.911 

7.938 

7.186 

6 !o02 

b.150 

5.799 

5.528 

5.319 

5.159 

42 

9.063 

8.065 

.292 

.692 

.227 

.866 

.586 

.370 

.205 

43 

,2l6j 

.191 

.398 

.782 

.304 

• 932 

.644 

.422 

.251 

44 

.368 

.318 

.505 

.873 

.381 

.999 

.703 

.473 

.297 

45 

.520 

.445 

.611 

.963 

.459 

6,066 

.761 

.525 

.343 

46 

9.673 

81571 

7.717 

7.053 

6753^ 

6.133 

5.819 

5.576 

5.389 

47 

.825 

.698 

.824 

.143 

.613 

.200 

.878 

.628 

.435 

48 

.978 

.825 

.930 

.233 

.690 

.266 

.936 

.679 

.481 

49 

.951 

8.036 

.323 

.767 

.333 

.994 

.731 

.527 

L5°_ 

9.078 

.143 

.414 

.845 

.400 

6.053 

.7  s 

.572 
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TABLE  3 (CONTINUED) 

orovisional  probits,  Y 
4.2  I 4.3  I 4.4  | 4.5  I 4.6 


5 • 1 37 
.175 
.212 
.250 
.287 


511  4.517 


5611  .560  .563 


807 

4.788 

4.778 

4.773 

4.773 

4.774 

837 

.817 

.805 

.799 

.798 

.799 

867 

.845 

.832 

.826 

.824 

.824 

897 

.873 

.859 

.852 

.8491 

.849 

.930 

4.913 

4.904 

■msEl 

4.900 

.959 

.941 

.931 

.926 

.925 

.987 

.968 

.957 

.952 

.950 

.015 

.995 

.983 

.977 

.975 

.044 

5.022 

5.009 

5.003 

5.000 
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TABLE  3 (CONTINUED) 


orovisiona]  probits,  Y 
- ! 5-3  I 5.4  I 5.5  j 5.6 


747  3.7401 3. 719 13.680 
772  tTt^T  TTpif  3.70S“ 
797  .790  .770  .732 

822  .816  .795  .758 
847  .841  .821  .785 


3-447  3.5o 
.675  .59 

.702  .62 

.729  .65 


.593  .482 

.621  .512 
.650  .542 


.42 2 3.272 
.452  3. 304" 
.482  .336 
.512  .368 
.542  400 


3.079  2.834 
3Til4~2TWl 
.148  .909 
.183  i .946 
.217  .984 


•iri; 


46 

4.900 

4.899 

4.895 

4.886 

4.869 

4.843 

47 

.925 

.924 

.921, 

.912 

.897 

.871 

48 

.950 

.949 

.946 

.938 

.924 

• 899 

49 

.975 

• 975 

.972 

.965 

.951 

.928 

50 

5.000| 5.000 

.997 

.991 

.978 

.956 

.562 

.600 

.637 

.675 

.713 
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TABLE  3 (CONTINUED) 


O'bS  . 

cfL 


\ % 

! p 

6.0  ! 6.1 

pi 

6.2 

-'ovisic 

6.3 

nal  probits, 

6.4  | 6.5 

Y 

6.6 

6.7 

6.8 

6.9 

0 

2.523 

2.132 

1.643 

1.030 

0.261 

1 

2 

3 

4 

2.564 

.606 
.647 
.688 
. 1 jyu 

'2.178 
. 224 
.270 

.316 

.362 

I.694 

.746 

.797 

.849 

.900 

1.088 

.146 

.205 

.263 

.321 

0.327 

.394 

.461 

.528 

.595 

7 

8 
9 

10 

2.771 

Oi  0 

4U±: 

.854 

.895 

.936 

2.408 

!i  r~  li 

.500 

.546 

.591 

1.952 

A /^A-^ 

£ . UVJ  P 

•055 

.106 

.158 

I.3B0 

). 

.*00 

.496 

.555 

.613 

0.661 
• ( dO 

.795 

.862 

.928 

0 067 

11 

12 

13 

14 

15 

2.978 

3.019 

.060 

.102 

.143 

2.637 

.683 

.729 

.775 

.821 

2.209 

.261 

.312 

.364 

.415 

1.671 

.730 

.788 

.846 

.905 

0.995 

1.062 

.129 

.196 

.262 

0.144 

.221 

.299 

.376 

.453 

16 

17 

18 

19 

20 

3 !l8^ 
.226 
.267 
.308 

.350 

2.867 

.913 

.959 

3.005 

.050 

2.467 

.518 

►-  /Q 
* ^ ^ 

.621 

.673 

1.963 

2.022 

.080 

.138 

.197 

1.329 

.396 

.463 

.530 

.596.. 

0.530 

.607 

.685 

.762 

.839 

21 

22 

25 

24 

2^ 

3.391 

.432 

.474 

.515 

.556 

3.096" 

.142 

.188 

.234 

.280 

2.724 

.776 

.827 

.879 

.930 

2.255 

.313 

.372 

.430 

.488 

1.663 

.730 

.797 

.864 

.930 

0.916 

.993 

1.071 

.148 

.225 

0.062 

.152 

.243 

.313. 

26 

27 

28 

29 

30 

3.598 

.639 

.680 

.721 

.763 

3.326 

.372 

.418 

.464 

.509 

2.982 

3.033 

.085 

.136 

.188 

2.547 

.605 

.663 

.722 

.780 

1.997 

2.064 

.131 

.197 

.264 

1.302 
• 379 
.457 
.534 
.611 

0.423 

.513 

.603 

.693 

.784 

31 

32 

33 

34 
36 

3.804 

.845 

.887 

.928 

.969 

3.555 

.601 

.647 

.693 

.739 

3.239 

.291 

.342 

.394 

.445 

2.83O 

.897 

.955 

3.014 

.072 

2.331 

.398 

.465 

.531 

.598 

1.688 

.766 

.843 

.920 

.997 

0.874 
• 964 
1.054 
.144 
.234 

0.050 

.156 

.262 

W 

37 

38 

39 

40 

4.011 

.052 

.093 

.135 

.176 

3-785 

.831 

.877 

.923 

.969 

3.497 

.548 

.600 

.651 

.703 

3.130 

.189 

.247 

.305 

.364 

2.665 

.732 

.799 

.865 

.932 

2.074 

.152 

.229 

.306 

.383 

1.324 

.415 

.505 

.595 

.685 

0.369 

.475 

.581 

.688 

■♦79-7. 

41 

42 

47 

44 

46 

4.217 

.259 

.300 

.341 

.383 

4.014 

.060 

.106 

.152 

.198 

3.754 

.806 

.857 

.909 

_^960 

3.422 
.480 
• 539 
■ 597 
.655 

2,999 

3.066 

.132 

.199 

.266 

27460 

.538 

.615 

.692 

.769 

1.775 

.865 

.955 

2.046 

.136 

0.900 

1.007 

.113 

.219 

..326_ 

0.035 

.162 

.289 

46 

4? 

48 

49 

50 

4.424 

.465 

.507 

.548 

.589 

4.244 

.290 

.336 

.382 

.428 

012 

.063 

.115 

.166 

.218 

3.714 

.772 

.830 

.889 

.947 

3.333 

.400 

.466 

.533 

.600 

2.846 

.924 

3.001 

.078 

.155 

2.226^ 

.316 

.406 

.496 

.586 

1 .432 
.538 
.645 
.751 
.857 

0.415 

.542 

.669 

.795 

.922 

1 
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